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A generalized linear sigma model for low energy QCD is employed to study the quark structure of 
eight low lying scalar isomultiplets as well as eight low lying pseudoscalar isomultiplets. The model, 
building on earlier work, assumes the possible mixing of quark anti-quark states with others made 
of two quarks and two antiquarks. No a priori assumption is made about the quark contents of the 
' states, which emerge as predictions. An amusing and contrasting pattern for the quark structure is 

I found; the lighter conventional pseudoscalars are, as expected, primarily of two quark type whereas 

, the lighter scalars have very large four quark admixtures. The new feature of the present paper 

^SJ ' compared to earlier ones in this series involves the somewhat subtle and complicated effects of SU(3) 

' flavor breaking. They do not alter the general pattern of two quark vs. four quark mixing obtained 

, in the SU(3)symmetric case but, of course, give a more detailed picture. 
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I. INTRODUCTION 

^\ 

Although it is an ancient subject, the study of light scalar mesons has greatly intensified in the last generation. 
^ I Some representative earlier works of this period are given in pj-p^. 

1^ I The "global" aspect of the scalar puzzle is the unusual spectroscopy of the light scalar nonet. At present, the scalars 
. below 1 GeV appear to fit into a nonet as: 



/ = : 77i[/o(600)] « 500 MeV 

>; 7=1/2: m[K] « 800 MeV 

/ = : to[/o(980)] w 980 MeV 

00 ■ 1=1: m[ao(980)] w 980 MeV (1) 

This level ordering is seen to be flipped compared to that of the standard vector meson nonet: 

^ ; 7=1: m[/5(776)] w 776 MeV nn 

O ■ 7 = 0: m[w(783)] w 783 MeV nn 

7=1/2: m[7i:*(892)] « 892 MeV ns 

7 = 0: to[0(1O2O)] « 1020 MeV ss (2) 

Here the standard quark content (n stands for a non-strange quark while s stands for a strange quark) is displayed 
at the end for each case. The vector mass ordering is seen to just correspond to the number of s-type quarks in each 
state. It was pointed out a long time ago in Ref. [l^], that the level order is automatically flipped when mesons 
are made of two quarks and two antiquarks instead of a single quark and antiquark. That argument was given for a 
diquark- anti diquark structure but is easily seen to also hold for a meson-meson "molecule" type structure which was 
advocated, at least for a partial nonet, in Ref. [2l|. Note that, in the "ideal" four quark picture, the states in Eq.((T]) 
consecutively have the quark contents: nnnn, nnns, nsns and nsns. One also notes that the masses of the putative 
scalar nonet are significantly lower than the other (tensor and two axial vector) p-wave quark- antiquark nonets. There 
are enough other scalar candidates [ao(1450), Ko(1430) and two of fo(1370), fo(1500,) fo(1710)] to make another nonet 
although the masses of its contents seem somewhat higher than an expected scalar p-wave nonet. Based on the usual 
effect that two mixing levels repel as well as some more detailed features, it was suggested [13] that a global picture 
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of these scalars might consist of a hghter "four quark" nonet mixing with a heavier "two quark" nonet. Further work 
in this direction has been presented by a number of authors [23j-(2a|. 

A field theoretic toy model to study these features was introduced in [2^ . Although simple in conception it turned 
out to be rather complicated to implement; as a result the present authors examined its features in further detail in 
(29I [sol . I31I [3^ . [33l . While the general treatment in ^29(1 only used the tree level Ward-type identities of the theory, the 
treatments in |3Cll | and (3lj | employed specific potential terms chosen in a systematic way to correspond to terms with 
the minimum number of underlying quark lines. Furthermore the quark mass terms were neglected, which not only 
is a simplifying feature but is also a check that the peculiar mixing results found are not an artifice of a particular 
choice of symmetry breaking terms. (In QCD it is expected that the main features of the various particle multiplets 
should also hold in the limit of zero light quark masses.) Next, in [s^l it was demonstrated that the results did not 
change much when the minimal mass term was included with equal masses for all three light quarks. It was also 
shown that the current algebra theorem (s^ l for pion pion scattering holds to a very good approximation in this more 
complicated case where the ordinary pion has a small admixture of a state containing two quarks and two antiquarks. 
The theorem has a very small correction due to a very small violation of the partially conserved axial vector current 
hypothesis in the model, however. Most recently, in [33] an amusing connection between the model and the instanton 
approach to QCD dynamics was discussed. This connection is not so surprising when one notices that it is the (broken 
by instantons) U(l)y!i symmetry which formally distinguishes the "four quark" from the "two quark" meson states. 

In the present paper we include the flavor SU(3) symmetry breaking in the model. As will be seen, this is not a very 
simple matter, even in the tree approximation being used. Within this framework we will not introduce any further 
approximations in our numerical treatment. Furthermore, the sensitivity of the predictions to the main uncertainties 
in the experimental inputs will be explicitly displayed. 

Considering the relatively large number of "outputs," it may be desirable to immediately just display our main 
"typical" results. These are the masses and the "two quark" vs. "four quark" percentages of the members of all 
four nonets (light and heavy pseudoscalars and light and heavy scalars) . They are listed in Tables [J and |TT1 Isospin 
but not SU(3) symmetry is being assumed. Note that for the 1=1/2 and 1=1 states, the prime denotes the heavier 
particle. For the 1=0 particles there are four states of each parity and they are denoted by subscripts 1, 2, 3, 4 in 
order of increasing mass. Altogether, considering the isospin degeneracy, there are 16 different masses. The 8 inputs 
comprise the pion decay constant, the four masses: [tott, ^tt', ma, "m'a ]i f^e strange to non strange quark mass ratio 
(which is related to assuming a value for rriK) and (as to be explained later) the sum and the product of all the four 
1=0 pseudoscalar squared masses for each of the possible scenarios for their identification with experimental states. 

Thus there are 9 mass predictions and 16 predictions for the two-quark and four quark percentages being displayed 
in Tables HI and [TTl To compare with Eq.lT]) for the "experimental" light scalar nonet, we read off the mass pattern: 

/i(742), «(1067), a(980) /2(1085), (3) 

in which only the a (980) mass was an input. Clearly, the pattern, featuring an approximate degeneracy of one 1=0 
state and the 1=1 state together with a flipping compared to the standard vector meson order, is very similar to the 
experimental one in Eq. ([T]) . Actually, these "tree level" predictions are expected to receive some non-trivial corrections 
as to be discussed later. We also read from Table HIl that the light scalar nonet is predicted to be predominantly of 
4-quark type. On the other hand, the heavy scalar nonet is predicted to be predominantly of 2-quark type and to 
have the standard (vector meson like) pattern: 

a'(1474), /3(1493), /«'(1624), /4(1784). (4) 

The conventional light pseudoscalar nonet is read from Table |T] as, 

^(137), X(515), r;i(553), 772(882), (5) 

where the structure is clearly close to the conventional one. Notice the well known fact that this nonet agrees in 
ordering with that of the standard vector nonet in the sense that the 1=1 state is lightest but disagrees in that the 
lighter 1=0 state is not at all close to the 1=1 state. At the theoretical level, this arises from axial U(l) (instanton 
type) terms as are included here. It can also be read off that the light pseudoscalar nonet members are predicted to 
be predominantly of "2-quark", i.e. qq, type. In contrast, the heavy pseudoscalar nonet members are predicted to be 
of predominantly "4-quark" type with a similar mass ordering as the light pseudoscalars; namely, 

7r'(1215), i^'(1195), r?3(1225), r?4(1784). (6) 

In this multiplet, the distortion of the pattern from that of the standard vector multiplet seems to reflect the greater 
role of "instanton" effects over "4-quark" effects. 
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A detailed discussion of how these resuhs were obtained is given in the foUowing sections. Section II briefly 
summarizes the "toy Lagrangian" and the choice of input parameters. The needed mass matrices are presented in 
section III and Appendix B. The first part of Section IV and, especially, Appendix A explain in a step by step way how 
the Lagrangian parameters (other than those associated with the "instanton" terms) are related to the experimental 
data. The second part of section IV explains predictions of the model and their sensitivity to changes in the input 
parameters. Section V explains the work associated with the instanton terms and the complicated 1=0 pseudoscalar 
sector. Finally, conclusions and some further discussion are presented in section VI. 
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TABLE I: Typical predicted properties of pseudoscalar states: qq percentage (2nd column), qqqq (3rd column) and masses 
(last column). 
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TABLE II: Typical predicted properties of scalar states: qq percentage (2nd column), qqqq (3rd column) and masses (last 
column). 



II. MODEL LAGRANGIAN AND PHYSICAL INPUTS 

The model employs the 3x3 matrix chiral nonet fields: 

M = S + i(t), M' = S' + i<j>'. (7) 

The matrices M and M' transform in the same way under chiral SU(3) x SU(3) transformations but may be dis- 
tinguished by their diff'erent U(1)a transformation properties. M describes the "bare" quark antiquark scalar and 
pseudoscalar nonet fields while M' describes "bare" scalar and pseudoscalar fields containing two quarks and two 
antiquarks. At the symmetry level with which we are working, it is unnecessary to further specify the four quark 
field configuration. The four quark field may, most generally, be imagined as some linear combination of a diquark- 
antidiquark and a "molecule" made of two quark-antiquark "atoms" . 
The general Lagrangian density which defines our model is 

jC = -^Tr (a^Ma^Mt) - ItV (a^M'S^M't) - Vo (M, M') - Vsb, (8) 
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where Vo{M, M') stands for a function made from SU(3)lx SU(3)r (but not necessarily U(1)a) invariants formed 
out of M and M' . 

As we previously discussed (30l |. the leading choice of terms corresponding to eight or fewer quark plus antiquark 
lines at each effective vertex reads: 

Vo = - C2 Tr(MMt) + Tr(MM^MM^) 

+ d2 Tr(M'M't) + e'^{eabce''"^MSM^^M'f + h.c.) 



C3 



detAf Tr(MM't) 



Tiln(:T-77+ ) + (1 - 7i)ln 



detAft^ ' ^ Tr(M'Mt) 



2 



(9) 



All the terms except the last two have been chosen to also possess the U(1)a invariance. The symmetry breaking 
term which models the QCD mass term takes the form: 

Vsb^~2Tt{AS) (10) 

where A = diag{Ai, A2, A^) are proportional to the three light quark masses. The model allows for two-quark 
condensates, aa = (5*°) as well as four-quark condensates Pa = {S'a)- Here we assume [sl] isotopic spin symmetry so 
Ai —A2 and: 

ai = "2 7^ as, = /32 7^ P3 (11) 

We also need the "minimum" conditions, 

m\./l>V.s\^„_ /S\ = 0. (12, 



dS / \ dS / ' \dS\ 

There are twelve parameters describing the Lagrangian and the vacuum. These include the six coupling constants 
given in Eq.Q, the two quark mass parameters, (^1 = ^2,^3) and the four vacuum parameters (ai = Q;2,a3,/3i = 
/32,/33).The four minimum equations reduce the number of needed input parameters to eight. 

Five of these eight are supplied by the following masses together with the pion decay constant: 

TO[ao(980)] = 984.7 ± 1.2 MeV 

TO[ao(1450)] = 1474±19MeV 

m[7r(1300)] = 1300 ±100 MeV 

= 137 MeV 

F„ = 131 MeV (13) 

Because TO[7r(1300)] has such a large uncertainty, we will, as previously, examine predictions depending on the choice 
of this mass within its experimental range. The sixth input will be taken as the light "quark mass ratio" A3/A1, 
which will be varied over an appropriate range. The remaining two inputs will be taken from the masses of the four 
(mixing) isoscalar, pseudoscalar mesons. This mixing is characterized by a 4 x 4 matrix M^. A practically convenient 
choice is to consider TrM^ and detM^ as the inputs. 

Given these inputs there are a very large number of predictions. At the level of the quadratic terms in the 
Lagrangian, we predict all the remaining masses and decay constants as well as the angles describing the mixing 
between each of (7r,7r'), {K,K'), (ao,ag), {k,k') multiplets and each of the 4x4 isosinglet mixing matrices (each 
formally described by six angles). 

Defining the total potential V = Vq + Vsb, the four minimum conditions explicitly read: 

d \^ \ 

=4e>i/33 + 4e;^/3ia3-2c2ai-h4c^a?-2Ai =0. (14) 



dSl 



dV 



8e;^/3iai - 2c2a3 -)-4c^a^ - 2^3 = 0. (15) 



dV 
dS'l 

dV 



Ae^aia3 + 2d2Pi =0. (16) 
= 4e^a2^2d2/33 = 0. (17) 
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III. MASS MATRICES 

In order to evaluate the eight independent parameters from experiment it is necessary to first obtain the formulas 

for the squared mass matrices corresponding to each set of particles with the same parity and isotopic spin quantum 
numbers. These are calculated by taking the second derivatives of the potential with respect to the appropriate fields. 
For the two mixing "pion" states (1=1, P=— ) we find: 



Mi 



4eg/33-2c2 + 4c2a? 4efa3 



(18) 



The 2 X 2 K meson matrix (1=1/2, P=— ) reads: 



(Ml) 



4eg/3i - 2c2 + 4cta? - Aclaia^ + Aclal Ae'^ai 
4 eg ai 2 d2 



(19) 



In the case of the two scalar a-mesons (1=1, P=+) the mass squared mixing matrix reads: 



-4e|/?3-2c2 + 12c4a? -4efa3 
-4 63 as 2 d2 



(20) 



The final 2x2 mass squared matrix describes the kappa-type scalars (1=1/2, P=+): 



-4 eg /?! - 2 C2 + 4 C4 + 4 ai 03 + 4 C4 al -4 63 ai 
-4egai 2d2 



(21) 



In the case of the 1=0 scalars there are four particles which mix with each other; the squared mass matrix then 
takes the form: 



4ef/35-2c2 + 12 c2af 
4\/2ef/3i 

4 el as 
4V2egai 



4x/2ei/3i 4eia3 4V2efai 
-2c2 + 12c2ai 4\/2e|ai 
4 eg ai 2 ^2 
2d2 



(22) 



For this matrix the basis states are consecutively, 



fa 

fb 

fc 

fd 



sl + sl 

%/2 

Si ss, 



^/2 



nn, 



nsns, 



(23) 



The non-strange (n) and strange (s) quark content for each basis state has been listed at the end of each line above. 
In the case of the four 1=0 pseudoscalars, the mixing matrix formula is a long one which is given in Appendix B. 



IV. PARAMETERS AND SOME PREDICTIONS 



A simplifying feature in the present model is that, as previously, the parameters other than C3 and 71 may be 
found without considering the four 1=0 pseudoscalars. This is due to the presence of "In's" in the last two terms of 
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the potential, Vq. So we consider the initial (six) parameters first. In Appendix A, a method is presented for their 
consecutive determination in a convenient way. 

Here we consider the strange to non-strange quark mass ratio, A3/A1 to be an input. For many years this has 
been known, from various analyses of the chiral treatment [36[, of the light pseudoscalars, to be around 25. In Fig[T] 
the dependence of the predicted decay constant ratio Fk/F^^ on the quark mass ratio, including the dependence on 
the parameter r7i[7r(1300)] is shown. For comparison, the experimental ratio [s^l is Fk/Ft^ k, 1.19. Thus, the model 
predictions are seen to be quite reasonable, although perhaps a bit too small for typical parameter choices. This ratio 
may be fine-tuned by including non-minimal kinetic terms in the Lagrangian. See e.g. Eq.(4.1) of ^38] . 

It is interesting to compare the values of the Lagrangian coefficients of the fully chiral invariant terms, 02,04,^2, ef 
as obtained here and displayed in Fig|5] with the corresponding ones obtained in the zero quark mass model: Fig. 2 



of [3l|. These parameters are clearly substantially similar to those in the zero quark mass case. This is a comforting 
feature since it agrees with the expectation that the light quark masses make only small changes in the QCD dynamics. 

On the lowest two rows of Figl2]we display the variations with A^/Ai and m[7r(1300)] of the two-quark condensates, 
the four-quark condensates and the quark mass type quantities Ai and ^3. 

In FigO we give the predictions for the masses of the strange mesons and the 1=0 scalars, showing the variations 
due to choosing different values of m[7r(1300)] and A^/Ai. It is seen that the mass of the ordinary kaon is better fit 
for A3/A1 closer to 30 than to 20. The predicted mass of the "excited" kaon is roughly similar to the (mentioned, 
but not established) candidate K(1460). 

The upper, right graph of Fig [3] displays the two predicted kappa-type particle masses. These are roughly consistent 
with the (discussed in [3I] but not established)i4rg (800) and the Kq{14:30). It should be mentioned that the masses 
obtained in this paper are being considered as "bare" ones, subject to non-trivial "renormalization" due to effects 
which provide unitarity corrections to the scattering amp litudes in which these particles appear as poles. This is 
illustrated for the single-M linear sigma model case in |19| . A similar remark applies to the particle widths obtained 
from the third derivatives of the potential. 

The graphs in the bottom row of Figl3] display the predicted masses of the four 1=0 scalars. Clearly the lighter two 
can be identified with the /o (600) [sigma] and the /o(980). Note that the very light sigma is an inevitable consequence 
of the present model and was not put in by hand. There are three candidates, /o(1370),/o(1500) and /o(1710)[one of 
which may be a glueball], for the heavier two particles. 

In FiglHthe predicted decay constants (i.e. the coefficients of the single particle terms in the appropriate axial vector 
or vector Noether currents) for the tt', the K, the K', the k and the n' are shown with their parameter dependences. 
Note that F-^i is very small but not exactly zero. This causes, as discussed in [sll, a very small violation of the 
"partially conserved" axial vector current ansatz in this kind of mixing model. 

One of our chief concerns here is the percentage "content" of "two quarks" (i.e. the qq content) vs. the percentage 
content of "four quarks" predicted to be in each particle state. In FigO these percentages are shown for the particles 
having non-zero isospin quantum number. For comparison with the SU(3) symmetric cases considered in [3l| and [3^ . 
note that the typical value of TO[7r(1300)] considered there was about 1215 MeV. Then we see that the ordinary pion 
has, as before, about a 15 percent 4-quark content (and 85 percent 2-quark content). The "excited" pion, 7r(1300) is 
predicted to have an 85 percent 4-quark content in this picture. Similarly, we see that the ordinary kaon is predicted to 
have about a 14 percent 4-quark content while the "excited" kaon has about an 86 percent 4-quark content. Notably, 
as before, the situation is drastically different for the scalar states. We read off that the lighter 1=1 scalar, ao(980) 
has about a 76 percent 4-quark content and the lighter kappa has about a 92 percent 4-quark content. 

The more complicated situation of the four mixing 1=0 scalars is described in FiglHl To read this, first note that 
the four physical states are labelled /i, /2, /a, A in order of increasing mass. Each is a linear combination of the basis 
states fa,fb,fc,fd given in Eq. (P5)l . fa and fb are 2-quark type while fc and fd are 4-quark type. Thus the lowest 
lying 1=0 scalar, /i = cr, taking j43/j4i=30 and m[7r(1300)] = 1215 MeV, has percentages in each basis state of, 

0.36, 0.04, 0.36, 0.24. (24) 

Altogether the sigma is about 40 percent 2-quark and 60 percent 4-quark in this case. This is similar to the roughly 
50-50 split we found in the SU(3) symmetric case (See Fig. 2 of^32]). /2, the next heaviest 1=0 state is about 95 percent 
4-quark and 5 percent 2-quark. The two heaviest 1=0 scalars both have majority two quark nature: /a is read to be 
about 63 percent 2-quark while is read to be about 93 percent 2-quark. 



V. 1=0 PSEUDOSCALARS 



The elements of the squared mass matrix for 1=0 pseudoscalars (given in Appendix B) are quite complicated and as 
a result it is convenient to treat this sector separately. As previously discussed, in this case the squared mass matrix 
contains two additional parameters (ca and 71) that only contribute to the properties of the 1=0 pseudoscalars. It is 
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therefore not possible to trade these two parameters with four experimental r] masses. Instead, we determine C3 and 
71 by fitting the trace and the determinant of the squared mass matrix to their corresponding experimental values, 
i.e. we solve for C3 and 71 from the following equations: 

[M^) = Tr (M^^)^^^ 
det (M^2) = det (m2)^^p (25) 

We identify the lighest two 77's predicted by our model (i.e. 771 and 772) with 77(547) and 77'(958) with experimental 
masses |36| : 

m'^'^P- [77(547)] = 547.853 ± 0.024 MeV, 

to'=''P [77'(958)] 957.66 ±0.24 MeV. (26) 

However, for the two heavier ?7's that our model predicts (i.e. 773 and 774) there are several experimental candidates 
below 2 GeV with masses [36j : 

m'^'^P- [77(1295)] = 1294 ±4 MeV, 
m'^'^P- [77(1405)] = 1409.8 ±2.5 MeV, 
w'^'^P- [77(1475)] = 1476 ±4 MeV, 

m'^^P- [77(1760)] = 1756 ±9 MeV. (27) 

We consider all six possible choices for identifying r/3 and 774 with two of the above four experimental candidates. This 
leads to six scenarios given in table Hill Equations (|25p result in a quadratic equation for 71 for which the discriminant 
versus 7Ti7r(1300) is plotted for all six scenarios in Fig. [T] We see that for A-^/Ai — 20 the scenarios 1, 2 and 4 are 
completely ruled out, whereas for A3/A1 = 30 up to to^(1300) « 1.25 GeV all six scenarios are possible. 

Therefore, for a given scenario there are two solutions for 71 and C3, and consequently two sets of predictions for 
the four 77 masses. We measure the goodness of each solution by the smallness of the following quantity: 

= r„-P(^fc) 



fc=i 



in which s corresponds to the scenario (i.e. s = 1 • • • 6) and / corresponds to the solution number (i.e. I = 1,2). The 
quantity Xsi x lOO gives the overall percent discrepancy between our theoretical prediction and experiment. For the 
six scenarios and the two solutions for each scenario, Xsi is plotted versus 77i7r(1300) in Fig. [H Clearly, scenario 3 is 
favored over the range of 77i7r(1300). 

Here, we present our predictions for the best fitting scenario 3 with A3/A1 = 30. (Of course, the other scenarios 
with different values of A3/A1 may also be of some interest for a more detailed look at additional properties of the 
system of the four etas). In the present case, the 77 masses are shown in Fig. [5] and we can clearly see a reasonable 
range of masses: The first mass is around 550 MeV (consistent with identifying it with 7^(547)), the second mass is in 
the range of 970 - 986 MeV (consistent with identifying it with 77(985)), the third mass is in the range of 1218 - 1250 
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FIG. 2: Parameters vs m[7r(1300)] (GeV): C2 (GeV^) (top left), (GeV^) (top right), C4 (second row left), 63 (GeV)(second 
row right), ai and as (GeV)(third row left), /3i and fis (GeV^) (third row right) and Ai and A3 (GeV^) (last row). (Note that 
d2, ai, /?i and Ai do not depend on the choice of ^.) 
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FIG. 3: Predicted masses (GeV) vs m[7r(1300)] (GeV): kaon system (top left), kappa system (top right), isoscalar salars (last 
row). 



MeV (consistent with identifying it with 77(1295)) and the fourth mass is around 1780 - 1790 MeV (consistent with 
identifying it with 77(1760)). 

The quark contents of the r;'s are given in Fig. [TO] and show a clear difference compared to those of the 1=0 scalars. 
The first two states (r;(547) and r;(985)) are dominantly of two quark nature whereas the two heavier states (r;(1295) 
and r;(1760)) are mainly four-quark states. 



Scenario 


T/'S 


1 


r;(547), 7?(958), 7?(1295), 77(1405) 


2 


77(547), 77(958), 77(1295), 77(1475) 


3 


77(547), 77(958), 77(1295), 77(1760) 


4 


77(547), 77(958), 77(1405), 77(1475) 


5 


77(547), 77(958), 77(1405), 77(1760) 


6 


77(547), 77(958), 77(1475), 77(1760) 



TABLE III: The six scenarios for identifying the predictions of our model for four 77s. 



VI. CONCLUSIONS AND DISCUSSION 

The results obtained [s^ seem to provide support for the picture of Hght scalars having a predominant content of 
two quarks and two antiquarks while the heavier scalars appear to be made mainly from one quark and one antiquark, 
as one would expect from p wave states in the non relativistic quark model. The predicted scalar and pseudoscalar 
masses were seen to be reasonably consistent with the experimental candidates listed in ^] . 
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FIG. 4: Predicted decay constants (GeV) vs m[7r(1300)] (GeV): pion system (top left), kaon (top right), kappa system (last 
row). The result for pion system is insensitive to Az/A\. 




While it appears a little unusual to think of say, the ordinary pion, as having some four quark content when treated 
in the effective Lagrangian framework, that is in fact the standard picture in the parton model approach to QCD. 
In the case of the two scalar nonets, the mass ordering itself naturally suggested such a picture. This picture was 
then inherited by the pseudoscalars when we chose to describe the scalars via a linear sigma model. There does not 
seem to be any problem with our treatment of the lighter pseudoscalar nonet. One might initially worry that the 
well established "current algebra" theorem for low energy pion pion scattering could get altered in the case of such 
a more complicated pion. However, we showed in [3l| and [s^ that there is only a small effect for this theorem. 
The situation concerning the heavier (mostly four quark) pseudoscalar nonet which appears in our model is neither 
so clear experimentally nor theoretically. Evidently this seems a fruitful direction for further investigation. The two 
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FIG. 6: Predicted quark contents of isoscalar scalar mesons vs m[7r(1300)] (GeV). The components [described in Eq. (|23|l ] are: 
fa (solid line), ft (dotted line), fc (dashed line) and fd (dotted-dashed line). 
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FIG. 8: Xsi vs m[7r(1300)] (GeV). 



new r/s which appear in this model are welcome from an experimental standpoint and the mixing with an expected 
glueball in this energy range is another interesting future topic. There is of course a similar glueball expected in the 
scalar isoscalar channel. 

In addition to the "global" analysis of the mass spectrum and substructure presented here, the "local" topic of 
the detailed structure of the pion and kaon scattering is also of great interest. Recently progress has been made in 
this area, strongly supporting the existence of both the sigma and kappa mesons, by using (^t'I and fi^l) the chiral 
perturbation theory approach combined with some dispersion theory results [4]| . Even though the present model is 
rather complicated, the fact that it is chiral symmetric means that we are guaranteed to get essentially the starting 
results of the chiral perturbation theory approach near the appropriate thresholds. In fact, since the light scalar fields 
are included, the starting results are even closer to the modern ones just mentioned. Thus we plan to next investigate 
the meson meson scattering scatterings in a number of different channels and at energies away from threshold. We 
are encouraged to go further by the results [l^ obtained some time ago by unitarizing the scattering amplitudes 
computed using a single chiral SU(3) nonet model. 

Actually, we have already started on this work and would like to briefly mention a result for the lowest lying scalar 
meson since it gives an idea of the accuracy to be expected for the masses we have presented in this paper. A more 
detailed analysis of the scattering amplitudes will be presented elsewhere. 

The mass spectrum of the light scalars receives considerable unitarity corrections due to low-energy rescattering 
effects. The masses obtained above appear as tree level quantities in the effective Lagrangian under discussion. 
Especially in the case of the scalars the physical states are rather broad and will appear as poles in the scattering of 
two pseudoscalar mesons. A simple way to estimate the scattering amplitude is to first compute the tree level scalar 
partial wave scattering amplitude and then unitarize it by using the K-matrix method. This is equivalent to an earlier 
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FIG. 10: Predicted quark contents of isoscalar pseudoscalar mesons vs m[7r(1300)] (GeV). The components [described in Eq. 
(|B11[) ] are: rja (solid line), rji, (dotted line), rjc (dashed line) and i]d (dotted-dashed line). 



approach [3| and amounts to replacing the tree level amplitude, Ttree by, 

^tree 



T = 



1 - iTf, 



(29) 



Following this procedure the sigma pole at 742 MeV discussed above appears in the unitarized pion scattering ampli- 
tude at z = AP ~ iMT with 



M = 477 MeV, V = 504 MeV. 



(30) 



This is of the same order as in [27|. Such scattering calculations should be performed to find the "actual" mass and 
width parameters for all the scalars in the present model. 
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APPENDIX A: PARAMETER EVALUATION 



By using the following formulas consecutively, it is possible to determine all the parameters of the model, one at a 
time, from the experimental inputs. The procedure is a more complicated version of the zero quark mass case given 
in Appendix B of [IH and the degenerate non-zero quark mass case given in the Appendix of jsl] . 

We first use the 1=1 scalar and pseudoscalar squared mass matrices to get the parameters. 



2d 



'2 



(«3e^f = ^iiml-ml,f-[^d,-{ml + ml,)r) (Al) 
The minimum equation then yields the ratio of the four quark to the two quark "condensate" : 



2(«3ei) 



(A2) 



Next, the tt — tt' mixing angle is found from the diagonalization of (M^): 

4:d2- ml - mL 



cos20Tr — 



^64 (aa e§)2 + 16 ^2 _ grf^ (ml + m^,) + {ml + ml,^ 
Then we get, from the formula for the pion decay constant (29j . 



(A3) 



1 

ai = - 



cos Ot, 



(^t) sine.' 



(A4) 



which then also gives /3i using Eq. (|A2p above. 

To go beyond this point we introduce the strange to non-strange quark mass ratio, A3/A1 as an input parameter. 
We may then solve for the quantity ai by using the two minimum equations: 

Al = 2P3ie3ai) + 2(3i{e^a3)-C2ai + 2clal 

A3 = 'iPi{e^ai)-C2a3 + 2clal (A5) 
Now, the rest of the parameters may be obtained by the consecutive use of the following equations: 



r'^ — 
C4 ~- 


1 2,2 
8 aj [ 


- ml 


- ml, 


- 16 


4C2 = 


16 alc1 + 4:d2- 


ml - 


ml, - 


ml - 


p"^ — 
63 — 


(eiai) 
ai 








as = 










/33 = 


ai Pi 
a3 









15 



{Ml) 



11 



-2d2 



I 2 : 2 2 2 \ 



8 



(Mi)l2 
(A^i)22 

COS 29 k = 



-4(e>i) 

2d2 



K)22 



(Mi)ll 



\/4(Afi)?2 + pfi)22-(A4)i 





(63 Q3) 

(e| Qi) 


j cos 9k 


_ ft 


1 1 (egai) 


sin 0/^- 




2 


cos 6*. — 


2 (It) Sine. 





(eg as)" (eg«3) 
(egai)2 (egai) 



(A6) 



With m[7r(1300)] = 1.215 GeV and A^/Ai = 30, these parameters are given in table HVl In this case, the rotation 
matrices that related the flavor bases to the physical states are given in (|A7p and 



C2(GeV'') 


1.62 


xlO" 


-1 


d2(GeV2) 


6.30 


xlO^ 


-1 


ef(GeV) 




1.68 






47.0 




ai(GeV) 


6.06 


xlO" 


-2 


as (GeV) 


7.68 


xlO" 


-2 


A (GeV) 


2.49 


xlO" 


-2 




1.96 


xlO" 


~2 




6.66 


xlO^ 


-4 


^3(GeV^) 


2.00 


xlO" 


'2 



TABLE IV: Calculated Lagrangian parameters:C2 , d2, ef, C4 and vacuum values: ai, Q3, /3i and /Ja, with m[7r(1300)] = 1.215 
GeV and A3/A1 = 30. 



0.923 0.385 
-0.385 0.923 



0.925 0.379 
-0.379 0.925 



0.493 
0.870 



-0.870 
0.493 



0.284 
0.959 



-0.959 
0.284 
(A7) 



0.601 0.199 0.600 0.489 

-0.107 0.189 0.643 -0.735 

0.790 -0.050 -0.391 -0.470 

0.062 -0.960 0.272 -0.019 



(A8) 



APPENDIX B: 

The elements of the symmetric matrix AI^ are given by: 

(M^) = (16 C4 ai/3i" + 16 C4 ai^as Pi P3 + 4: C4 ai^'aa^/Ja" - 16 63 ai^^i^^s - 16 63 ai^aa Pi P^'^ 

-4 63 ai^a^'^P^^ - 8 C2 ai^/^i" - 8 C2 ai^ag Pi P3 - 2 C2 ai^ag'/Jg' - 16 C3 7i'ai'/3i' 
-32 C3 7i^ai as /3i /^s - 16 C3 Ti^as^/Js" - 32 C3 71 ai^/3i" - 32 C3 71 ai as /^i /Js 

-16 C3 ai^pi^) I ((2 ai Pi + as /3s)' ai^) (Bl) 
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(m2) ^2 = (-4 ^P^ (4 63 ai^as + 4 eg ai^as'/^i^/^s + e3 ai a3^/3i + 4 c3 7i^ai2/Ji2+ 

4 C3 7i^Q!i Q!3 /?! /^s + 4 c3 71 OLx'^x + 2 C3 71 ai as /Ji /^s + 2 cs 71 as^/Js^ + 2 C3 ai as P1P3)) 
/((2ai/3i+a3/33f aias) (B2) 

(M^) = (-16 63 ai^as - 16 63 ai a3^/?i /?3 - 4 63 a3^/?3^ - 16 C3 7i^ai A - 16 C3 71^03 /?3 

+16 C3 71 as /33 + 16 C3 ai / ((2 ai Pi + as Psf^ (B3) 

K)i4 = (-4V2(4e3ai^ + 4 63 ai^a3 Pi Pz + 63 a-i^as^Ps^ + 2 C3 7i^ai as /3i 
+2 C3 7i^a3^/93 - 2 C3 71 a3^/33 - 2 C3 ai az P\)) 

/{2aiPi+azPzfai (B4) 

iK)22 = (l6c4aiW/3i^ + 16c4aia3^/?i/?3 + 4c4a3^/?3^ -8c2ai^a3^/?i^ -8c2aia3^/3i/33 
-2 C2 az^pz'^ - 32 C3 7i^ai^/3i^ - 32 C3 71 "3 P1P3-8 cz az^Pz^) 

/ ({2aipi+azPzfaz'') (B5) 

(^^)23 = ("'^ ^ ("^ ^3 ai^as /3i^ + 4 63 ai az^Pi Pz + 63 03=^ /Sg^ + 4 C3 7i2ai /3i - 4 C3 71 ai Pi 

+2 C3 71 as /33 - 2 C3 az Pz) ^i) / ((2 ai Pi + az Pzf (B6) 

iK)2A = (^1 - 1) (2 71 "1 A + as Ps) cz) I (2 ai pi + as /3s)' (B7) 

(M2)gg = (8 ai V^i' + 8 rf2 ai az Pi Pz + 2 ^2 ag'/^s' - 16 C3 71'ai' + 32 C3 71 ^i' - 16 C3 ai^) 

/(2ai/3i+a3/33f (B8) 

(^')34 = (^1 - 1)' ai ^3) / (2 ai pi + az Pz)'' (B9) 

(M^)^^ = (8 d2 ai'^pi^ + 8 d2 ai as PiPz + 2 c«2 as'/3s' - 8 cs 7i'q!3' + 16 cs 71 az^ - 8 cs as') 

/i2aipi+azPz)' (BIO) 
The basis states for the above matrix are: 



V2 ' 



rib = 03> 



V2 ' 

% = <^3'- (Bll) 
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C3(GeV^) 


-3.78 X 10"* 


71 


5.27 xlO~^ 



TABLE V: Calculated parameters: cs and 71. 

With m[7r(1300)] = 1.215 GeV and A3/A1 — 30, these paprameters are given in table fVl 
In this case, the rotation matrix for 1=0 pseudoscalars becomes: 

" -0.675 0.661 -0.205 0.255 ' 

_ 0.722 0.512 -0.363 0.291 

° -0.134 -0.546 -0.519 0.644 

0.073 0.051 0.746 0.660 
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